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I. INTRODUCTION 


The linear instability of small amplitude disturbances in boundary- 
layer flows has received considerable attention in an attempt to develop a 
better understanding of this aspect of the transition process. This 
instability is sufficiently described by the parallel stability theory, 
where the boundary-1 aver flow is assumed to be parallel. In spite of the 
qualitative success of this assumption, experiments have shown systematic 
differences from the parallel stability theory. 

Apart from predicting a critical Reynolds number that is lower than 
that given by the parallel stability theory, evidence from the experir.;ents 
show that the growth rate of the disturbance is a function of the coordinate 
normal to the wall in a subsonic or supersonic boundary layer. The parallel 
stability theory cannot predict that and consequently a detailed comparison 
with experiment cannot be made in a meaningful way to the necessary degree 
of precision, the parallel theory gives only a qualitative results. Also in 
cases where large pressure gradient or large amount of suction exists as in 
laminar flow control systems, the parallel theory fails to predict 
accurately the stability characteristics. 

With the advancement of the numerical techniques and the possibility of 
solving the disturbance equations to considerable accuracy, a more critical 
comparison with experiments seems now possible. This demands a consistent 
theory describing the behavior of the disturbances over large regions of a 
real growing boundary layer. The nonparallel stability theory, which takes 
into account the dependence of the flow parameters oh the chordvise 
coordinate (see Figure 1 for the coordinate system), as well as the 



velocity normal to the the v/all, this theory gives a more accurate 
mathematical model for the development of the disturbance. The full non- 
parallel incompressible stability problem for tv/o-dimensional boundary-layer 
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flows has been solved by Bouthier *, Nayfeh et al , and Raster , and 

5 - 8,22 

applied to different types of boundary-layer flows . The nonparallel com- 
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pressible problem has been reported by El-Mady , El-!-Iady and Nayfeh 

11,12 13 

for two-dimensional flows, and by El-Hady and Nayfeh for three- 

14 

dimensional boundary layer flows. Tlie work of Gaponov , which appeared 
during writing this report, investigates only the nonparallel effects on 
two-dimensi onal di sturbances. 

It is the purpose of the present work to describe a formally correct 
method, based on the nonparallel linear theory, that examines the two- and 
three-dimensional stability of compressible boundary- 1 ayer flov/s. The 
method is applied to the supersonic flat plate boundary layer at Mach number 
4.5, and the theoretical growth rates of disturbances are compared with the 
experimental results of Kendall The method is applied also to the 
i nf i ni te- span swept wing boundary-1 ayer flow with suction at Mach number 
D.82 to evaluate the effect of the nonparallel flow on the development of 
cross-flow instabilities. 
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II. NOrJPARALLEL SPATIAL STABILITY THEORY 


This study is concerned with the stability of small amplitude 

disturbances in two types of boundary -layer flows, the flat plate supersonic 

boundary layer, and the infinite- span swept wing transonic boundary layer 

with suction. In both types, the mean flow is independent of the spanwise ^- 

coordinate. However, the most ar.iplified Tol 1 mi en- Schl i cht i ng instabilities 

can be chordwise (two-dimensional) or oblique (three-dimensional) for the 
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flat plate supersonic boundary layer , but certainly oblique for the 
infinite-span wing transonic boundary layer. In the latter, crossflow 
instabilities also exist due to the presence of boundary-layer crossflov/. 

He consider disturbances that is generated by a steady source 
oscillating at the frequency co* (monochromatic wave) at the curve x* = 0. 

This spatial wavetrain is assumed to be of uniform amplitude in the 
spanwise coordinate, but with variable phase. With this in mind, we present 
here a formulation for the stability problem of the infinite-span wing 
boundary layer. All spanwise variation of the meanflow and disturbance 
quantities is zero. 

We begin with the equations of motion for compressible flow, and 
introduce dimensionless quantities by using as reference, the values at the 
edge of the boundary layer, and a reference length L* defined as 

L* = (v4x*/U^)l/2 (1) 

where u* is the kinematic viscosity coefficient, x* is the distance 
perpendicular to the leading edge, and U* is the velocity in the chordwise 
X* - direction. Here, * denotes a dimensional quantity, and e denotes 



the condition at the edge of the boundary layer. Then, the Reynolds number 
is defined as 

R = U*L*/ue = (UexVM^/^ ( 2 ) 

We consider the boundary -layer flow to be weakly nonparallel, that is 
the meanflow quantities are weakly varying functions of the chordwise 

coordinate, and the normal velocity component V is small compared to the 
velocity components U (chordwise) and W (s panwise). The meanflow variation 
is scaled by introducing the new variable X defined as 

X = ex , e = 1/R 

then we write the mean-flow velocity field as 

U(X,y), V(X,y), W(X,y) 

and the mean-flow temperature, density, viscosity, and pressure fields as 

0 (X,y), p(X,y), y(e), P(X) (4b) 

We express the mean-flow quantities in terms of mean and disturbance 
motions. The disturbance motions are assumed to vary slowly with the 
chordwise coordinate. They are described by an amplitude which is uniform 
in the spanwise coordinate and a phase function <})(x, 2 ,t) that vary with 
both x and z coordinates as well as time t. We use the method of multiple 
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scales and expand the disturbance components in power series of e. Then 
the mean-flow quantities in dimensionless form can be written as 


( 3 ) 

! 

\ 

\ 

(4a) 
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u = U(X,y) + [uo(X.y) + f:Uj(x,y)] exp (5a) 

V = V(X,y) + Cvo(X.y) + eVi(X,y)] exp (i4>) (5b) 

w = W(X,y) + [wo(x,y) + eWi(X,y)] exp (i(|)) (5c) 

P = P(X) + [P„(X,y) + ePi(X.y)] exp (U) • (5d) 

p = p(X,y) + !>q (X,y) +epi(X,y)] exp (i<J> ) (5e) 

0 = e(X,y) + [6o(X,y) +ce^(x.y)] exp (i«l> ) (5f) 

p = y (o) + “-[e^(X,y) + (x,y)] exp (i4> ) (5g) 

Where the real frequency is 

to (X) = - (t>t (6) 

and the complex wavenumber components are 

»(x) = *x (n 

e(x) = *2 (8) 


Here the subscripts t, x, and z denote partial differentiation with respect 
to these variables. Equations (7) and (8) define a complex wavenumber 
vector as 

K(X) = grad 4> (9) 

and hence 

Curl K = 0 (10) 

We substitute Eqs. (3)-(8) into the dimensionless equations of 
motion, and transform the spatial and temporal derivatives from x, z, 
and t to X snd 4 > . We eliminate nonlinear terms of the product of the 
disturbance quantities since they are assumed to be small, subtract the 
average equations of the mean-flow, and equate terms order c" and order c. 
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The eigenvalue problem that is formed by the zero-order equations can 


be written as 


[a] {Zq>= 0 (11) 

with the boundary conditions 


Zoi = Zo3 = Zo5 = Zo7 

o 

It 

at y = 0 

(12) 

Zoi » Zo3 , Zo5, Zo7 

0 

as y -»■ " 

(13) 


where 

Zoi= U o» Zq 2= LI ojr* Zq3 = = P, 


^05~ ®o* 


"06 


= 6 


oy 


07 


= w, Z 


08 


= W 


oy 


(14) 


and the elements of the matrix a are defined in Appendix I. The boundary 
condition (13) assumes that the disturbances are subsonic. It is reminded 
here that stability formulation is presented for the infinite-span-wing 
boundary layer, where Eqs. (11) consist of eight first order system of 
differential equations. The same set of equations governs the three- 
dimensional stability of the flat plate boundary layer, except that the 
spanwise mean velocity component W will disappear from the elements of the 
coefficient matrix a of Eqs. (11). Hov/ever, the two-dimensional stability 
of the flat plate boundary layer is governed by a set of six order systems 
of differential equations. 
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Equations (11) - (13) are integrated nuirerically to provide a solution in 
the form 

{Zq}= A(X) {^(X.y) } (15) 

where c is an eigenrnode solution vector that is normalized in some specified 
but arbitrary manner, and A is a complex amplitude function which is yet to 
be determined. 

Using Fq. (15), we write the non-ho::iogeneous first-order equations as 
{ Z } - [a] { = ib^)A +{b 2 )A 

J A 


' ^ 3 * 2,5 ” = 0 at y = 0 


Zjj, ® ^sy 


(17) 

(18) 


where the elements of the vectors bj and b 2 are given in Appendix II. They 

are known functions of a, 3, C, and the meanflow quantitites. By using the 

* 

adjoint function C(X,y), which is a solution of the adjoint homogeneous 
problem; 


(t*y } + [a] ic*} = 0 


= c* = t* = = 

0 

at 

II 

o 

4 6 8 




• C* ?*, ^ 

4 6 8 

0 

as 

y 


we write the solvability condition of Eqs. (16) - (18) as 
Qb^} A +(b^}A^] Ic*/dy = 0 


(19) 

( 20 ) 
( 21 ) 


( 22 ) 
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where T indicates the transpose of the vector. We can write Eq. (22) in the 
following form for the evolution of A: 


Ax = i a(X)A 


(23) 


where 


” * T /” * T 


(24) 


The solution of Eq. (24) can be written as 


A = A^exp [ij «(X)dX] = A exp [ie J“(X)dx] 


(25) 


where Aq is a constant. 

The evolution of the wavenunber components aandgare detemvined as follows. 
We replace Z with Z in Eqs. (11) - (13), differentiate the result with 
respect to X and use the symmetry condition (10) to obtain 


- Ca] } = i {b2} + {d } 


(26) 


with homogeneous boundary conditions similar to those given in Eqs. (12) and 
(13), the vector d is defined in Appendix III. Applying the solvability 
condition to Eqs. (26) we obtain 




oo 

= - J{d } U*)'^ dy 


0 



0 


{ b2}{C*}^dy 


(27) 
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Similary, differentiating Eqs. (11) - (13) with respect to Z and using the 
syinmetry condition (10) we obtain 

= 0 (28) 

The evolution of “ and 3 is governed by Eqs. (27) and (28). Equation (28) 
indicates that the complex spanwise wavenumber 6 will be constant along the 
chordwise coordinate. 

Tiierefore, to the first approximation, Eqs. (5), (17), (25), and (27) 

gi ves 

{ Z 0 A 0 ^ J'(ot + ) dx + i 6z - iojt] (29) 

The effect of the growing boundary layer is reflected in Eq. (29) by 
perturbing a to gS and by making the eigen solution function of the 

chordwise coordinate in addition to being function of the normal 
coordinate. In Eq. (29) both Band to are constants. 
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III. NUMERICAL PROCEDURE 


One can write the real and imaginary parts of a and B as function of a 
real wave number vector with direction y , and a real spatial growth rate 
vector a with direction that is 


Re(“) = k cos '1' , Im(“) = -ocos? 


(30) 


Re(B) = k sin'l' , Im(6) = -osin'l' 


( 31 ) 


For a given R, we specify Re(B) (or the angle f) , Im(B),and a dimensionless 
frequency defined as 


F = 2 Tif u* 


(32) 


v;here f* is the dimensional frequency in CPS, and “ denotes freestream 
conditions. In all numerical results presented in this paper, we chose 

Im (e) = 0 at the leading edge, that is = o. We note that Eq. (28) 
requires that im(B) will remain constant along the chordwise coordinate. 

The eigenvalue problem (11) - (13) is integrated numerically from y = 

ye to the wall using a variable step size algorithm based on the Runge- 

Kutta-Fehlburg fifth-order formulas, coupled with an orthonormalization 
1 8 

procedure. Following the same procedure, the adjoi nt problem (19)-(21) is 
solved for the adjoint fuction c*, except that it has. the same eigenvalues 
as the homogeneous problem (11) - (13). 

To determine S'(x), vve need to evaluate a^^d “x* “x 

calculated using Eq. (27), while evalijated by solving the homogeneous 

part of Eq. (26) with homogeneous boundary conditions. 

ll 



Using Eq. (29), the spatial growth rate of the disturbance, defined as 
Re(Z^/Z), can be written as 

a = -lin(a) - e [Im( a ) - Re ( 33 ) 

The first term is the spatial growth rate in a parallel flow, while the second 
is a correction due to nonparallel effects. The dependence of the eigen func- 
tion C in the chordwise coordinate in the nonparallel theory, made it possible 
to determine the growth rate as function of the normal coordinate as well as the 
disturbance flow variable. 

Because a hot wire anemometer in a supersonic stream responds to both 
mass-flow and total temperature disturbances, we use these disturbance flow 
variables to demonstrate the theoretical results. The amplitude function of the 
mass-flow disturbance is given by 

pu = (34) 

and the amplitude function of the total temperature disturbance is given by 

T = C, + (y- 1) m! U (35) 


12 



IV. FLAT PLATE BOUNDARY LAYER AT M =4.5 

oo 

In this section, we compare the theoretical growth rates of a 
three-dimensional disturbances with Y = 55°, and a two-dimensional distur- 
bance ( V = 0°) in the flat plate supersonic boundary layer with the 
experimental findings of Kendal 1^^. 

A. Theoretical Results 

For a threo-diirionsional disturbance with wave angle 4'= 55° and frequency 
6 

F = 30 X ICT (first mode). Fig. 2 shows theoretical results of the varia- 
tion of the mass-flow amplitude |pu| as well as the total temperature 
amplitude ] T | with normal coordinate y at R = 1550. Both are normalized 
with IpuI the maximum value of |pu^ Tlae variation of |puI and |t| with y 
is typical for the unstable frequency range at R = 1550, except that the 
y-location of the peak value varies slightly with frequency. The peak loca- 
tion of jpu 1 > for example, changes from y = 12. 7 at F = 5 x 10"® to y = 12.0 
at F = 75 x 10"^, the peak value is not sharp in general, but extends over a 
considerable normal distance as shown in figure 3. While | T ] has two peak 
values, ]pul has only one peak that is twice as much as the higher peak 
value of lx) and occurs at different y-locations. 

The normal structure of the different amplitude functions evolves 
differently as the disturbance travels downstream. Hence, at each normal 
location, lpu| and |T| vary differently in the chordwise direction. Because 
of this behaviour, growth rates based on the amplitudes I Pul and It I are 
functions of y at fixed x-locations. This is shown in Fig. 2 as o1pu| 
and o| T | respecti vely. The discontinuity in o|pul curve at y = 6 is due to 
a near zero amplitude of 1pu| at this y-location. 
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Figure 2, Variation of the amplitudes and growth 
rates of the mass flow and total temperature 
disturbances with the normal coordinate at 
R = 1550 (first mode). 




I I I I I I 
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Figure 3. Normal location of the peak value of the 
mass-flow amplitude at R = 1550 and different 
frequencies (first mode). 


15 



At the same chorcVise location, R = 1550, we calculated the growth 
rates of a two-dimensional disturbance at higher frequencies (second mode). 
For the frequency F = 130 x lOr ^ , fig. 4 shows the variation of the mass- 
flow amplitude jpujand the total temperature amplitude I. tI (normalized v/ith 
the maximum value of IpuI ), together with the variation of the growth 
rates o|pu| ando] T | with y-coordi nate. 4gain, the variation of |pu| and 1 T | 
with y is typical for the unstable frequency range at R = 1550, with the 
y-location of the peak value varies slightly with frequency. The peak 
location of |pu| varies from y = 12.40 at F = 120 x 10"^ to y = 12.05 at 

-6 , I 

F = 158 X 10 . Figure 5 shows the variation of the peak location of lpu| 

with y for this frequency range, and shows that wide peaks exist at some 
frequencies. In constrast with the first mode, fig. 4 shows that for the 
second mode the peak value of I T I (the higher one) is twice as much as the 
peak value of |pu| and occurs at different y-locations. 

B. Comparison With Experiment 

It is worth noting that, like the theory, the experiment also gives 
different stability conditions depending on the normal location of the 
measurements as well as the disturbance variable measured. 

Here we compare the theoretical results with the experimental findings of 
Kendall This experiment possesses some advantages that made it the most 
suitable for comparison with the linear stability theory. The experiment 
was performed in the JPL 20-inch wind tunnel where the tunnel wall boundary 
layer were laminar. A small disturbance of almost pure frequency was 
produced by a glow-discharge generator, and introduced at an angle 0° and 
f = 55“ . The chordwise growth rate of the disturbances were determined for 
boundary -layer flow at Mach number 4.5 by following the maximum amplitude 
point. Three-dimensional disturbances, introduced at f = 55“, were found to 
grow nearly equal at various spanwise z-locations. 
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Figure 5. Normal location of the peak value of the 
mass-flow amplitude at R = 1550 and different 
frequencies (second mode). 
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A hot wire in a supersonic flow responds to both mass-flow and total 
temperature disturbances with a ratio depending on the sensitivity coeffici- 
ents that are determined by mean-flow conditions and the mean temperature of 
the hot wire. Unless measurements are conducted following a contour of 
constant y-location (y*/L* constant), sensitivity of the hot wire to changes 
in the mean-flow conditions cannot be avoided. Also, the hot wire response 
when operated at high constant overheat is proportional to the mass-flow 
disturbance, its sensitivity to the total temperature disturbance is 
mi nimized. 

In Kendall experiment, the distribution of the measured disturbance 
aiiplitudes across the boundary layer were almost similar at different 
x-locations. Hence, the maximum amplitude points followed in the experiment 
had a contour of almost constant y-location, and sensitivity of the hot wire 
to changes in the mean-flow conditions was minimum. The hot wire was 
operated at high constant overheat to minimize its sensitivity to total 
temperature disturbance. Hence, energy fluctuations obtained by traversing 
the hot wire through the boundary layer at different x-locations for 
particular frequency, have been interpreted as being nearly proportional to 
the mass-flow disturbance amplitude. 

Figure 6 shows theoretical growth rates of a three-dimensional 
disturbance with wave angle f = 55° (first mode) as function of the 
dimensionless frequency F at R = 1550. Parallel results, that is neither 
function of the normal coordinate nor the disturbance flow variable, is 
indicated by a dashed line, it coinsides with the results of Mack given in 
reference 19. Nonparallel growth rates, indicated by a solid line, are 
calculated applying Fq (33) at the y-location where the mass-flow 
disturbance has a maximum. Kendall measurements of the growth rates 
following the maximum of a disturbance amplitude detected by the hot wire 
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are given also in the same figure for coirparison. Nonparallel results 

enhances the comparison between the theory and experiment by shifting the 

theoretical ly predicted growth rates to higher values compared to parallel 

results. However, higher growth rates are given by the experiment in the 

_6 

frequency range from 10 - 35 x 10 

Figure 7 shows the theoretical growth rates of a second mode 
disturbance with = 0°, as function of frequency at R = 1550. Parallel 
results indicated by the dashed line compares with that of Mack given in 
reference 19. Nonparallel results, indicated by the solid line, are calcu- 
lated applying Eq. (33) at the y-location v/here the mass-flow disturbance 
has a maximum. Experimental growth rates by Kendall are given also in the 
same figure for comparison. The maximum of the disturbance amplitude 
detected by the hot wire is followed in these measurements. Nonparallel 
results again enhances the comparison between the theory and experiment by 
shifting the theoretically predicted unstable frequency range to higher 
values compared to parallel unstable frequency range. 

It should be pointed out that the location of the peak detected by the 

hot wire at fixed x-location was assumed to be the same for all frequencies 

2 0 

as indicated later by Kendall , snd that slightly different experimental 
data could have been attained for other frequencies by searching out the 
peak corresponding to a particular frequency as we did in the theory. The 
spread of the experimental data shown in Fig. 6 and Fig. 7, and the diffi- 
culty in comparing it with the theoretical results is expected. This 
experiment, like others, is not designed specifically to verify the non- 
parallel stability theory. Hence, the effect of the nonparallel flow, prob- 
ably, has contributed to the difficulty in comparison. The dependence of 
the hot wire response not only on the mass-flow disturbance but also on the 
temperature disturbance, even with a minimum unknown percentage, has 



Figure 7. Variation of the growth 
two-dimensional disturbance, 

mode) with frequency, at R = 


rates of a 
0° (second 
1550. 
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contributed also to this difficulty, figure 2 and figure 4 show that the 
variation of the growth rates ^Ipul and ^ Ix I are different with the 
y-location. 

The wide peaks of |p u | at some frequencies make it difficult 
experimentally to measure the actual growth rates, specially if growth rates 
change rapidly around the peak location as indicated in Fig. 2 and Fig. 4. 

This is demonstrated in Fig. 8 which shows the variation of the growth rate of 
the mass flow disturbance with the normal coordinate around its peak amplitude 
at the frequencies F = 15 x 10“^ (first node), and F = 125 x 10"^ (second 
mode). Experimental growth rate measurements around these frequencies are 
sho\cn in Fig. 8 to fall within the wide peak of the amplitude |pu | that is 
off 1 .5 -3 percent from the peak value. Changes in the growth rates of I P ul 
corresponding to an average of 2 percent deviation from the peak value is 
calculated for all unstable frequencies of the first and second modes. 
Calculated growth rates form a band that is shown by the hashed area in 
Figs. 6 and 7. Most of the experimental data fall within this band. 
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Figure 8, Variation of the growth rate of the mass 
flow disturbance with the normal coordinate 
around its peak amplitude for a first and a 
second mode frequencies, and the corresponding 
experimental growth rate measurements. 
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V. BOUNDARY LAYER ON A TRANSONIC SWEPT WING 


The nonparallel stability theory presented in section 2 is applied to the 
boundary layer with suction on a 23° swept infinite span wing. The airfoil sec- 
tion is supercritical with a normal chord c = 1.98 m, and the freestream Kach 

number is 0.82. Parallel stability characteristics of the boundary layer on 

IP 21 

this wing have been investigated by El-Hady and Mack . However, seme incom- 

12 

plete nonpai'allel calculations have been reported by El-Hady for the same wing. 

Figure 9 shows the distribution of the pressure coefficient Cp on the upper 
surface of the sicfoil section together with the distribution of tne suction coefiicient 
along the chord. The suction coefficient is defined as = -PcVc./p«U,o, where the 
subscript ^ denotes wall condition. 


On a sweptback wing the boundary layer is three-dimensional due to the 
presence of boundary layer crossflow. Because the crossflow profiles are highly 
dynamically unstable, crossflow instabilities generate and dominate both the 
leading edge and rear parts of the wing where the mean crossflow components have 
large values. 

We advance a band of initial wavelengths of constant frequency disturbances at 
the leading edge. These disturbances develop downstream according to Eqs. (25), 

(27), and (28) . The spanwise component of the wavelength ^z/c (normalized with the 
normal chord) remains constant, and 7 remains zero if its initial value at the 

leading edge is zero. We calculate the parallel and nonparellel growth rates as well 
as the logarithmic amplitude ratios of these wave components as they develop downstream. 


The logarithmic amplitude ratio N is calculated from 


where R^ is the initial 
stability curve, and a^ 


° ( 36 ) 

Reynolds number at the lov/er branch of the neutral 
is the corresponding disturbance amplitude. The o- i n 
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the above fornula is calculated by using Eq. (33) for a constant frequency F. 

All results presented here are for the stationary crossflow disturbances, that 
is F = 0. 

Calculations at the leading edge region show that a wave component with 
constant ^z/c =0.00054 has the largest anplitude ratio, however, a wave com- 
ponent with constant ^z/c = 0.0036 has the largest amplitude ratio at the rear 
part of the wing. For these wave components, fi gure 10 shows the variation of 
IpuI and o)pu| with y at the chordwise locations x/c = 0.015 (for ^z/c = 
0.00054) and x/c = 0.90 (for ^-z/c = 0.0036). The growth rate alpu| changes 
rapidly around the peak, of 1pu|. We use the value of a|pu| at the y-location 
where |pu| is maximum to represent the nonparallel results as we did in the pre- 
vious section. 

Figures 11 and 12 show the parallel and nonparallel unstable regions of 
three wave conponents around the most unstable one (that gives the largest 
amplitude ratio) as well as their growth rates as they develop downstream. At 
the leading edge region, figure ll i ndicates that the nonparallel maximum growth 
rates are 30 - 45 percent higher than the parallel. It is worth noting that 
this region extends to about 10 percent of the chord. The boundary layer is 
accelerated rapidly from the leading edge to x/c » 0.025, and small suction 
rates are applied to the boundary layer starting at this location. It is clear 
from figurellthat in this region, parallel and nonparallel maximum growth rates 
are affected only by the rapidly accelerated boundary layer rather than suction 
rates. 

At the rear part of the wing, figure 12 shows that the nonparallel maximum 
growth rates are 15 - 30 percent less than the parallel. The boundary layer in 
this region is severely decelerated and high suction rates are applied to it 
(peak of Cs 0.0018). The sudden drop in the parallel growth rates around 
x/c = .84 is due to the peak suction rate. Here, maximum growth rates are 
affected by both the adverse pressure gradient and the severe suction rates. 
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Figure 10. Variation of the amplitude and growth 
rates of the mass flow disturbance with the 
normal coordinate for the most unstable wave 
components at the leading edge and aft regions 
of the wing. 





I 



Figure 11. Growth rates of most unstable wave 

components of constant leading edge 

region, as they develop downstream. 
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Figure 12. Growth rates of most unstable wave 
components of constant X^/c at the aft region 

j t A 


as they develop downstream. 




Figure 13 shoivs the variation of the maxintjm logarithmic amplitude ratio 
v/ith the wave components of constant >z/c* At the leading edge, 

the boundary layer flow is destablized due to the nonparallel flow 
effects. Higher nonparellel is shown with nearly 64 percent maximum dif- 
ference compared to the parallel contrary, at the rear part the 

boundary layer is stabilized due to nonparallel flow effects with a maximum 
reduction in of nearly 43 percent. Figure 11 shows also a slight shift in 

Ilia X 

the value of the most unstable wave component ^z/c due to nonparallel flow 
effects. Kie '’"'ost unstable wave component at the leading edge reaches its maxi- 
mum amplitude at x/c » 0.032 (R = 900) while at the rear part, the maximum 
amplitude is reached at the trailing edge. 
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'max. 



Figure 13. Variation of the maximum logarithmic 
amplitude ratio with wave component of constant 
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VI. CONCLUDING REMARKS 


1. A method is presented to investigate the two- and three-dimensional 
stability of compressible boundary layer flows. The method is suitable for the 
flat plate boundary layer as well as the boundary layer on an infinite swept- 
back wing, where the meanflow is independent of the spanwise coordinate. 

2 . No'nparallel stability results for the flat plate boundary layer at Mach 
number 4.5 are in better agreement with Kendall experiments than the parallel 
stability results. 

3. Nonparallel flow effects on crossflow instabilities are calculated for 
a specific infinite span transonic swept wing. The destabilizing effect at the 
leading edge and the stabilizing effect at the rear part is believed to be 
directly affected by the pressure gradients and the suction rates at this 

regi on. 
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APPENDIX I. 


The nonzero elements of the matrix a in equation (13) 

ai2 = 1 
^21 " 

^22 ^ 


P 


^23 

^24 

^25 

^26 


= RU ' - ia y ' [jn + ■) J 
P0 M 

= iaR (m + 1 } vM^a$ 

p 

= Cm + 1 ) ^ - (fiU ' ) ■ 

0 V 


831 = -lot 


^33 ^ — 


0 

^34 ' 


®35 ” — 

0 


0 


®37 


341 = -ixa C2p' ^ ro') 
u 0 

®42 ” 

843 = x( rp'0' ^ r^ _ i^ . - 6^) 

P0 0 ~ p$ " 

2 

®44 = “■'Aryl'l^ CeI + ^) 4 ’ + ciU' + ew 



P 0 

(i+ r) (aU' + ew) + rp'$ 
P 0 pG 


are 
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^46 " 

0 

947 = -1x3 r^) 

P 0 

^48 " ■■'^® 

®56 " ^ 

852 ■= - 2 r Cy - I) H%' 

853 = Rre' . 2tr (y - 1) CaU' + eu'l 

]iQ 

- i3L iy - ^ 

P 

arc = ^ ^ ^^'co 15 . CU- + W- I - h!!. 

-p 

355 = z^ilL 

p 

agg - -2r (y - 1 1 nfw< 

®78 " ^ 

^83 p1 + Cm + 1) ^ + RJ;!^ 

V 0 ">10 

®84 “ iHg _ Cni + 1 ) ByM^3> 

p 

^85 ~ ^ ^ ^ CPW ^ ) * 

0 ]3 

as6 “ 

p 

ap7 - + “^ ■*■ 

P0 

388 " 

p 

where 


$ = all + &W - u 

X = 1/CR + ir yM^ 4 >) 

p 


37 



r; > 


= ^ 
do 

r is Prandtl number, m = 2 Ce - l)/3, and r = 2 (e + 2)/3, 
corresponds to the Stokes hypothesis, it is taken e = 0.8 i 
The prime denotes differentiation with respect to y, and p 
for the meanflow viscosity. 


where e= 0 
n this analysis, 
is used here 
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APPENDIX II. 


The nonzero elements of the vectors and in equation (18) are: 


bi2 = C^I^ + Cgig 

^13 " S^l 

^14 ^ ^4^1 ^6^3 

^15 " ^7^5 

*^18 " S^l ^2^4 


and 


^22 " 


^2^7 

^23 " 

^ 3^6 


^24 " 

^4^6" 

^5^6 ^ CgIg 

^26 = 

S'lO 


II 

CO 

CM 


C 2 I 9 


where 

= -I Cm + 1) a©, C2 = -R, C3 = 0, C4 = rx Cu'0 + 20), Cg ^ rx0. 


C5 = Bx> ^-7 = -Rr. Cg = -i (m + 1 ) 30 

V y 

and 



- vMf 

0 


Uy + V - 2 CU0y + V0') 
^ 0 ^ 


Uy + V - 1 CU0y + V0') 
^ 0 ^ 


5 + UCgx + VCg + 1^0^ 
0^ 

54 - 2 ^^ t‘^‘ 4 X - 

0 


*=1 

h 

0 ^1X 
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I« = ir (ay)y - 1 Uy C, + ( 2 ir ya - U) +1^ (m + 1 ) 
2 R~ ^0^ '“R“ g'^R 


(C3x + ^'3 ^ 


3X 


yVll (UUj^ + VU') ?4 - V Cg - c 


0 


0 


4X 


+ ia{i (rUy + mV) + 1 (UUy + VU') - i* 3 y Wy ?c + im Puy ?? 
~R~ ^ ^ 

1-3 “ i Coty)y “ y 5 o (2i yoi ~ U ?oy Cm 1 ) p ^?y 

+ yx ?2 + (m + 1) yU'j^ + U'y^ + ryV" Cg + yU ' Cg + (rV + mU^) X 
Cy'?g ■*■ P ?g) my ' ?ix ~ y. ^ 3 


I 4 = (1 gPx - i '’1 ^ i Cm + 1) p 3 ^IX + (ay)x Cj 

R 0 ^ 

- Cuw^ + VW) CyM^ ■ ‘’5 ■ '^X '*’ -— ■*■ 


+ (21 pa - U^) ?7y “ y. ?o 
R e e ^ 


I5 = Cy - 1 ) + 2tp Cac^ + 6C7) CrU^ + mV + W^) - i X 

R 0 

©X Cl + 2 Cy - 1 ) U' + CmU^ + rV ) V 3 - CU0x 

R 

+ ve') Cih! C4 - tj ) * (y - 1) m2 * ''''4> ^ ^ K * 


« 7 

'5 * ‘If '5X - I '6 

'e ■ ■! Cc, + ym|u C4 - U 4 g) 

I, = (2ir ycx - U) ;, + 1 y' + (m + 1) p (Vo + ieCj) -C 4 

■y -T e ' W 

Iq = 1 mp'c, + Cm + 1 ) y + yU' Cc + Cn pa - U) Co 
^ R ‘ ' ° R 0 
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1 


In = i(m + 1) + ( 2i VO - y.) 

^ R ' R 0 

1,0 = (v - 1) li (Ml C 3 + U 5 ^) + (|i VO - U) 
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APPENDIX III. 


The nonzero elements of the vector d in equation (32) are: 

U'R (T_ - B^) + ia (B, + m B^) 


IT 


<^2 “ ^1 ^2 
V0 

- (iaR + mayM^B^) + maBg - 1_ (pU" + 1 CpU')^ 

" >' p2 


+ PX - ®6 > '6 


d-5 ~ Bp + iB^ C 


2 '=4 '"5 ■’5 


d^ =» ‘‘“^x’^7 ^7 ^^®3 *^1 ''' ^“^x^7 ^2 " ^ Ctp' + 

■p 


r0 * B, - tBq + Bo + :xD-| B^) Co + i>"XYMo 

0 ^ P0 


°4 - ^2 '^2 + 


(Bj + B^)» - B,p + lx Dg B,d Oj r + D3 

' '■p X 030 


■ m CBo + ®vl - X Bp D 

0 -P0 ^ 


'7 *^6 


?5 + trx (Bg - |i B7) ?5 + iBx X 


2 


C2B3 + rB^ + Bj) Cy + TBx By Cg 


dg = ”^8^*X ^2 ~ (Br_ B^ + Rr0 p^ + 2iDg B-|q) c^ '*’ ^ Bg p^ - B 2 ) C^ 


P“0 


+ Bg - 2 p Dg (U'U'j^ + W’W'^) - Dg (U’^ + W'^) ,pv + B^, 

P0 p ^p^X 


■'' ^®3 ^6 “ ^ ^6 '^'x ^8 


dg = iB (B3 + m B4) + ^ (W'j^ - B^W) C3 - (IM B2)X 

P0 

54 + mg (B2 - i ©x) - i CCW'j^)' - 1 (W'p^)' - 1 p„ CaW' )'- C5 

0 0 p p ~~2 

V 

- Cw^ P.)x Cg + iR Bg Cj + B3 Cg 
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where 


h " i 1 ®X 
y 0 

^3 " ly ^X ■ 1 ^'x 

p 

B. = 0' Oy - 1 0 'y 

^ ^ ^ 0 

85 ' fz - ®X 
0 0 ^ 
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y 
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D, = 2y ' + re* 
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D„ = r (y - 1) M‘ 
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